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We extend the analytical arguments of Metlitski and Grover (arXiv:1112.5166) to compute the 
entanglement spectrum and entanglement entropy of coplanar antiferromagnets with SO( 3) order 
parameter symmetry. The low-energy states in the entanglement spectrum exhibit the tower of 
states structure, as is expected for systems that undergo spontaneous symmetry breaking. Our 
results are consistent with numerical results on the triangular and Kagome lattice. 


PACS numbers: 05.30.-d,75.30.-m 

A system that exhibits spontaneous continuous sym¬ 
metry breaking (SSB) in the thermodynamic limit N —> 
oo, where the order parameter does not commute with 
the Hamiltonian, will for any finite system size have a 
unique ground state. For example, a finite size quantum 
antiferromagnet will have a completely symmetric spin 
singlet as a ground state. 

How then can any finite system show signatures of its 
symmetry breaking fate? Anderson 1,2 pointed out that 
such a finite size system contains a ’tower of states’ or 
’thin spectrum’ 3,4 : eigenstates with an energy 0(1/N) 
that vanishes in the thermodynamic limit, and who have 
a degeneracy structure that reveals the order parameter 
symmetry. This tower of states can be used to discover 
the symmetry broken state numerically 5 ’ 6 , which deals, 
by its very nature, with finite systems only. 

Recently, the notion that the ground state entangle¬ 
ment can reveal the low-energy spectrum of a phase has 
gained widespread recognition.' For SSB-systems, this 
suggests that information about the symmetry breaking 
lies hidden in the completely symmetric ground state. 
Indeed, it has been shown numerically for the collinear 
antiferromagnet that the entanglement entropy does ex¬ 
hibit the tower of states by having an extra logarithmic 
term. These results were confirmed analytically using 
the non-linear sigma model. 9 

In this paper we extend the analytical results of Ref. [9] 
to explain the numerical results 10,11 for the entanglement 
structure of coplanar antiferromagnets on the triangular 
and Kagome lattice. Unlike collinear antiferromagnets, 
the order parameter now has an SO(3) symmetry, which 
changes the degeneracies of the tower of states. In Sec. 

I we will introduce the concept of SO( 3) order, and the 
corresponding tower of states. We will then show in Sec. 

II that, consistent with the numerical results, the tower of 
states will also be present in the entanglement spectrum 
of the ground state. Our results support the thesis that 
symmetry broken systems can be understood by studying 
their finite size unique symmetric ground state. 


I. COPLANAR ANTIFERROMAGNETS 

In this section we review the main relevant properties 
of coplanar antiferromagnets, which will lead up to the 
representation of coplanar systems in terms of an S'O^) 
nonlinear sigma model. 

A. Triangular antiferromagnet 

As a prime example of a coplanar antiferromagnet, let 
us consider the Heisenberg model 

h = jY j s 1 -s j , ( 1 ) 

(ij) 

with J > 0 and the sum over nearest neighbor sites, 
on a triangular lattice. On a triangular lattice, the unit 
vectors that connect neighboring sites are 

e 1 = (1,0), 

= < 4 ^' 

e* = (2) 

Unlike its square lattice companion, the triangular an¬ 
tiferromagnet is frustrated already at the classical level. 
Nonetheless, a classical ground state exists where the tri¬ 
angular lattice is split into three sublattices A, B and C, 
and spins on different sublattices make a 120 degree an¬ 
gle which each other. For example, we choose an explicit 
reference state (see Fig. 1), 

n A = (1,0,0), 

= (-^,- 7 ,^, 0 ), 

nc = (-^\/3,0), (3) 

such that our reference ordered state will have (S.^a) = 
an A , etcetera, where a is the spin of a single site. 

The order parameter is an element of the rotation 
group SO( 3): each ordered state can be obtained from 
the reference state of Eqn. (3) by a unique rotation 
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FIG. 1: (color online) The triangular lattice with the classical 
120 degree coplanar ordered state. This configuration is used 
as reference state, see Eqn. (3). 

R £ 50(3) acting on the reference state. This makes 
the triangular antiferromagnet coplanar , in contrast with 
the square lattice antiferromagnet which is collinear and 
where the order parameter is just a vector lying on the 
2-sphere. 

The question is whether the ground state of the quan¬ 
tum triangular antiferromagnet bears any resemblance to 
the classical 120 degree state. As early as 1973 Ander¬ 
son proposed that the ground state of Eqn. (1) is not 
ordered, but instead a resonating valence bond ’spin liq¬ 
uid’ state. 12 Contradicting numerical results (for example 
Ref. [13] versus Ref. [6]) keep this an open problem, but 
for the present purpose, we assume that the model does 
exhibit symmetry breaking. To avoid the subtleties asso¬ 
ciated with a possible spin liquid phase, one can stabilize 
the 120-degree ordered state by adding a small ferromag¬ 
netic coupling between spins on the same sublattice. 

Similarly, coplanar order can be found in Kagome an- 
tiferromagnets. 

B. Tower of States 

On any finite size system, however, the ground state of 
the triangular antiferromagnet has not broken the spin 
rotation symmetry. After all, the total spin 5 = Y • Si 
commutes with the Hamiltonian and thus the ground 
state is an eigenstate of total spin too. It turns out the 
ground state is given by the symmetric singlet (S = 0) 
state. 27 

Insights into this singlet ground state and its 
low energy excitations can be gained from separat¬ 
ing the Heisenberg Hamiltonian into a spin-wave part 
and terms that depend on the total spin on each 
sublattice. 1-6 To separate the short-wavelength from 
long-wavelength physics, we perform a Fourier transform 

o 1 &ikri q 

_ y/N e 

fl = j./^ 7 A'S- k (4) 

k 

where 7 k is a geometric factor associated with the lat¬ 


tice: 7k = Y <5 e * fc<5 where <5 are the four six neighbors 
on the triangular lattice as defined in Eqn. (2). At the 
momenta k = (0, 0), k = K = (27r, 2n/y/3) and k = —K 
the geometric factor | 7 k| reaches the maximum value of 
1 which prevents a spin wave treatment at these momen¬ 
tum points. Isolating this part of the Hamiltonian, 

H 0 = 3 J (sLo - 5 K ■ 5_ K ) (5) 

and expressing it in terms of the total spin S = Y; Si and 
the spin on each of the sublattices, Sa = Y^ieA Si, etc., 
we uncover the triangular equivalent of the Lieb-Mattis 14 
Hamiltonian, 

Ho = “(Sj-Sb+Sb’Sc + S.i’Sc) 

= % {&-%-%-%)■ <«> 

This Hamiltonian can be solved exactly. 5,6 Its ground 
state is such that the spin on each sublattice is maximal 
and equal to oN/3 where N is the total number of sites 
of the system, and the total spin of the system is zero: 
the aforementioned singlet state. The excited states have 
a higher value of the total spin, 

9 7 

E=— 5(5+1) +E 0 . (7) 

These states are referred to as the ’tower of states’ or the 
’thin spectrum’. 1-6 The lowest energy excitation is to the 
total triplet state, with energy 0{1/N ) above the ground 
state. The eigenstates with energy Eqn. (7) are (25'+1) 2 
degenerate, since the states consist of combinations of the 
three sublattice spins. 5,6 

In the thermodynamic limit, the tower of states be¬ 
comes degenerate with the singlet ground state and the 
symmetry broken state can be formed as a superposition 
of the tower states. 

The remaining part of the Heisenberg Hamiltonian 
contains short-wavelength physics at momenta k ^ 0, K. 
It contains linearly dispersing spin waves following Gold- 
stone’s theorem. Since to = c|k| for small but nonzero 
momenta, and the lowest possible momentum on a finite 
N = L x L size system is fc min = , the lowest state 

with an excited Goldstone mode in d = 2 has an en¬ 
ergy 0(l/\/N ) above the ground state. Therefore the 
tower states are energetically separated from the spin- 
wave states. 

This separation allows for numerical detection of the 
tower of states, which was done by Bernu et al. 6 for lat¬ 
tices up to 36 sites big. They indeed found that the low 
energy states displayed the structure of Eqn. (6) - thus 
enabling to see the precursor states of spontaneous sym¬ 
metry breaking in a finite size system. 

C. Nonlinear sigma model 

The Heisenberg antiferromagnet is amenable to a semi- 
classical treatment, in which the Hamiltonian is mapped 
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onto the action of the so-called nonlinear sigma model 
using spin coherent states. An introduction to this map¬ 
ping for square lattice antiferromagnets is given in Refs. 
[15] and [16]. 

For the triangular antiferronragnet, the mapping has 
been first derived by Dornbre and Read 1 '. Let us quickly 
summarize their result. The spins on each site are 
parametrized by coherent states |JV), with the property 
that they are eigenstates of S\ with eigenvalue the vector 
N. In imaginary time, the action becomes 


[f* 

S = ia dr 

Jo J 

[ du^Ni- 

o 

(dNi 

x 

dNA 
dr J 


rR 

+a 2 J 

Jo ' 

(a) 



(8) 

For the anticipated 120-degree 
parametrize these vectors as 

coplanar 

order, 

we 

Ndr) = -= 

R(t) (di + aLi(r)^j 


(9) 


y 1 + 2aLi(r) ■ Hi + a 2 L 2 (T) 

where Hi are the reference vectors Eqn. (3), that only 
depend on the sublattice. The rotation matrix i?,(r) is 
an element of 50(3), and Li{r) expresses the average 
magnetization, and a is the distance between neighboring 
sites. 

We then expand the action for small L , take the con¬ 
tinuum limit and integrate out the average magnetization 
fluctuations of L. An action that depends solely on the 
coarse-grained 50(3)-field R(x,t) emerges, 


q _ rs 

“ "T. 

+Tr 


d 2 x I dr ( —Tr 


(R~ L d T R)~ 


P{R~ l VRy 


( 10 ) 


where P a p = § b c n ia n ifi is the projection opera¬ 

tor onto the plane of the magnetic order of the reference 
state. This projection causes the spin waves to have dif¬ 
ferent velocities depending on whether their polarization 
is in the plane of the order or out-of-plane. 

The nearest-neighbor triangular Heisenberg model 
gives rise to the 50(3) non-linear sigma model of 
Eqn. (10). In general, based on symmetry argments, 15 
an 50(3) model can also contain terms of the form 

Let us now dis- 


Tr 

P(i?- 1 a r i ?) 2 

and Tr 

(ir^vi ?) 2 

cuss these terms : 

in full generality. 


D. Charts on SO(3) 

It seems that the model Eqn. (10) depends on the ref¬ 
erence state, which should not be the case for a symmet¬ 
ric theory. However, by choosing the appropriate chart 
on 50(3) we can reveal the correct anisotropy structure. 


The group of rotations in three dimensions has three 
degrees of freedom, which can be parametrized in vari¬ 
ous ways. Starting with the standard basis for the corre¬ 
sponding Lie algebra so(3) of infinitesimal rotations, 

/0 0 0 \ / 0 0 1 \ /0 -1 0 \ 

T 1 = 0 0 -1 , T 2 = 0 0 0 , T 3 = 1 0 0 , 

\0 1 0/ V-l 0 0/ \° 0 0/ 

(n) 

we can parametrize elements of 50(3) using the expo¬ 
nential map 


R = R(Co) = e a ' f (12) 

where |w| < 7r. The derivatives R~ l d ll R in the model 
Eqn. (10) can be expanded for small rotations |w| <C 1, 

Tr {R-'dpRf = -2 48 ^ , 2 M ^cc) 2 

M 

“ jj|2 - 4sin2 

~ —2 (d^u) 2 + ... (13) 

where in the last line we have expanded for small u> up 
to fourth order. Notice that small |w| <C 1 corresponds 
to small rotations. Similarly, the anisotropic part can 
be expanded as Tr P(R~ 1 d fJ R)' 2 « —(c^c o x ) 2 — (d^Uy) 2 — 
2 {d^ z ) 2 . 

A more insightful parametrization is through the rota¬ 
tion quaternions: four real numbers a,, i = 0,».. 3, with 
the constraint a 2 = 1. A rotation matrix as a function of 
quaternions is 

/ o 3 \ 3 

Rij = Sij I a/ - '^2 Ofcflfc ) + 2 aiOj + 2 ^2 £ijkaoa k (14) 
V k =l / k -l 

and the integration measure is the same as for the 3- 
sphere, 


lSO( 3) 


dR = 


1 

27r 2 



(15) 


The isotropic term in the nonlinear sigma model now 
becomes 


Tr(R~ 1 d fl R) 2 = -8 (d^a) 2 . (16) 

The anisotropic term becomes quite elegant in this nota¬ 
tion. The anisotropy direction should, of course, depend 
on the average magnetization, but the model action has 
to be symmetric. This can be done by introducing a 
smooth vector field on the 3-sphere, v(a) = that 
defines the local anisotropy direction. With the specific 
choice of vector field direction, 


( 0 0 0 1 \ 
0 0 1 0 I 
0 -10 0 
-1 0 0 0 / 


(17) 
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we can use this anisotropy vector field to write 
Tr(l - PXiT^i ?) 2 = 4(a 7 1 c) M a) 2 . 


(18) 


Therefore, the most general SO( 3) nonlinear sigma 
model action is 


S =\ J d ^ x J q dT {x\\(d T a) 2 +pi\(Vd) 2 + (x±-X\\)(a'r 1 d T a)‘ 2 + (p ± -p ll )(a'y 1 Va)‘ 


(19) 


where \ is the spin susceptibility and p is the spin stiff¬ 
ness. Renormalization group studies 18-20 argue that any 
anisotropy (that is p± py or \± ^ X||) is danger¬ 
ously irrelevant close to the quantum phase transition to 
the disordered phase. Furthermore, a Holstein-Primakoff 
(large S) expansion has been performed around the 
coplanar ordered state , 21-23 finding anisotropy in both 
the susceptibility as well as in the spin stiffness. 

Notice that the quaternion chart is equivalent to the 
chart of two complex numbers Zi with \z-if 2 = 1 to repre¬ 
sent rotations . 15,19 


II. ENTANGLEMENT SPECTRUM 

Having introduced coplanar antiferromagnetism and 
the most general action that describes it, we proceed by 
analytically computing the entanglement spectrum of the 
ground state of Eqn. (19). We thereby follow the ap¬ 
proach of Ref. [9], which uses the nonlinear sigma model 
to derive the entanglement spectrum of a collinear mag¬ 
net. 

The entanglement entropy has been computed numer¬ 
ically by Kolley et al. for the triangular lattice 11 as well 
as for the Kagome lattice 10 . We find that our analytical 
results are consistent with their numerical results. 


A. Ground state wavefunction 

First we need to find the ground state wavefunction 
of the model Eqn. (19). Therefore we split the k = 0 
component from the spin wave fluctuations, 

a(x,r) = n 0 (r ) \/l - TTi(x,T) 2 /pi + e.^fx, r)/y/pi (20) 


and so the action becomes, up to quadratic order, 

S = J dT(d T n 0 ) 2 

+ \ J J dT (^(^ T7r *) 2 + (^ r *) 2 ) ( 22 ) 

where c 2 = Pi/x-i represents the spin wave velocity 
squared. 

We will now transform back to the Hamiltonian for¬ 
malism. The Hamiltonian associated with the zero- 
momentum part no(r) is that of a free particle on a 3- 
sphere, 


H t = 


Lip 

2 X||V 


(23) 


The ground state has zero angular momentum, which is 
a singlet, and has a wavefunction that is constant on the 
3-sphere. 

The Hamiltonian for the spin-wave part, that is k ^ 0, 
is 


H - = £ < 24 > 

k^0,i 

where the dispersion et — C|||fc| for the inplane modes 
i = 1,2, and e| = cj_|fc| for the out-of-plane mode. Bear 

in mind that Cj = \J~Pifxi- The static propagator that 
should be reproduced by the ground state wavefunction 
is 

j(x)) = D i:j (x , y) = Sij^ ^2 ( 25 ) 

fc^O * 


where e) with i = 1,2,3 are unit vectors orthogonal to 
?fo(r). The anisotropy is parametrized by pi = p\\ for 
i = 1,2 and p± for i = 3. Notice that the 3-direction 
is defined by the anisotropy vector field 7 1 no ■ In or¬ 
der to keep the same number of degrees of freedom, 
the 7 r,; fields cannot have zero momentum, which means 
J d 2 xni(x,T) = 0 . The anisotropic term expanded for 
small 7 Ti yields an extra term for the 3rd spin wave field, 

( 21 ) 


and hence, with the inverse of D defined as 

J d 2 zQij (x , z)D jk (z, y) = ^d 2 (x - y) - 5 ik (26) 

we find the ground state wavefunction of the spin-wave 
part is 

i/j[iri] oc exp f d 2 xd 2 yTTi(x)Qij(x,y)-!Tj(y)\ (27) 


{aXd^a) 2 = (< 9 ^ 3 (: t , t )) 2 + ... 
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The contribution from spin-waves is anisotropic, depend¬ 
ing on the direction of the magnetic order, characterised 
by no G 50(3). One must thus bear in mind that the 
degrees of freedom 7 r,(a?) are associated with a specific 
choice of the ordering direction. In particular, if we 
choose the global rotation no = (1,0,0,0), the corre¬ 
sponding 50(3) element is the identity, and hence the 
order the same as in the reference state: in the spin xy- 
plane. In this case, we can identify the 3rd spin wave 
held n 3 (x) with the actual z-component of the spin. 

For further convenience, we introduce the inverse cor¬ 
relator Q(x,y ) without dependence on the spin-wave ve¬ 
locity, 

Q{x,y) = Y j 2\k\e i ^ s -V. (28) 

Uo 

and the spin wave anisotropy matrix 
/ c ll 0 0\ 

c = 0 C|| 0 = cyl + (c_l — cy)-P 3 . (29) 

\0 0 c_l/ 

where P 3 projects onto the third axis. With these defini¬ 
tions, the ground state wavefunction for the 50(3) held 
is 

ip[no,-Ki\ oc exp J d 2 xd 2 yKi(x)cijQ(x, . 

(30) 

B. Reduced Density Matrix 

Our next step is to compute the reduced density ma¬ 
trix pa- The density matrix of the total system is dehned 
as p(a,a') = tp[a\^*[a']. In order to obtain the reduced 
density matrix, the total system is split into region A 
and region B, with volume Va and Vb, respectively. We 
will give the degrees of freedom on each region a corre¬ 
sponding subscript, so a(x) = cla{x) + Sb(x). In order 


to get the reduced density matrix on region A, that is 
Pa = Trsp, we must identify the two helds on region 
P, that is as{x) = a' B (x). Because the ground state 
wavefunction is Gaussian and thus deviations from the 
average magnetic order are exponentially suppressed, we 
can expand a and a' both around a common choice of 
directions, and subsequently integrating out ttb- Later, 
we need to restore the broken symmetry of choosing a 
direction around which to expand. 

Let us expand the helds around the north pole of the 
3-sphere, 

a°(x) = \Jl- TTi(x) 2 /pi (31) 

a\x) = TTi(x)/\J~Pi- (32) 

The ground state wavefunction is still Eqn. (30), but now 
the helds 7 r-i(x) are unconstrained and can have a zero 
momentum component. In fact, the magnetic order of a 
held configuration n l (x) is given by the k = 0 component, 

no = e 0 + y J d 2 XTr l (x)ei/^pi (33) 

up to hrst order in 7 t 1 (x), where eo = (1, 0, 0, 0), and e) • 
e 3 = Sij form a four-dimensional orthogonal basis. Notice 
that since / d 2 xQ(x, y) = 0, the inclusion of the zeroeth 
momentum component does not change the ground state 
wavefunction. 

It does allow us, however, to relate 7 f l B (x) = 7 t l B \x). 
We introduce the short-hand notation 

kaQab'kb= / d 2 x d 2 yit 1 A Qij(x,y)Tt J B , (34) 

J A J B 

where Qbb is the matrix operator where the hrst and 
second coordinate are both in region B, and so forth 
for similar expressions. The reduced density matrix now 
becomes 


p A oc s Vi tb exp 
1 


aQAA'X A + 'x’aQaA'x’a + (aA + 'x'a^Q AB^B + 'XbQBa(k A + 7T ’a) + 2'KbQBB^B^J 


oc exp 


^27t aQaa~xa + 27t 'aQaa'x'a ~ ( n A + ^'a)QabQbbQba{^a + 'x’a) 


(35) 


where f B d 2 zQ BB (x, z)Qbb{z, y) = 5 2 [x—y) 1. Notice that the anisotropy matrix-structure of Qaa and QabQb^Qba 
are the same, namely that of the spin-wave matrix anisotropy matrix Eqn. (29). Using the scalar definition of Q(x, v) 
of Eqn. (28), it follows that 


Pa[tta^a] oc exp 


— — ( 2 'kacQaa'ka + 277(4 cQaa'x'a — ( 7r a + t^'a)^ abQ bbQ ba^ a 


tt'a)) 


(36) 
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We can write this result in terms of the sum and difference of tta and tt' a 
Pa[tta,tt' a \ oc exp 


— O ( 7r ^4 ~ Tt'a)cQ AA {tT A — 7 T 'a ) — o ( 7r ^4 + 7 r ^l)c ( QaA ~ QabQbbQba) Qt A + TT'a ) 


(37) 


r 


The anisotropy in this result is relative to the expansion 
point we chose in Eqns. (31)-(32). In order to restore 
the symmetry we split off the k = 0 component of the 
774 -fields from the k 0 components. That is achieved 
by a Fourier transform the fields and Q-functions. The 
difference term in the reduced density matrix reads 


with the usual commutation relation [n a , pp\ = i5 a p. 
On the 3-sphere, the total angular momentum squared 
^2a</3 ^a /3 * s e d ua l to the Laplacian p 2 so that for large 
I we have 


(n|e 21 | n') « (n\e 21 \n') 


exp 


1 


-o vl J2 ( wa ~ K' A y(k)c ij Q{k)(TtA - 7T' A y{-k ) 


( 38 ) 

Recall that in order to derive this expression, we as¬ 
sumed that tta and tt' a are close to each other, so that 
the local anisotropy vector field can be approximated by 
+ n'). The anisotropic part of the wavefunction, 
which defines the orthogonal direction, can thus be writ¬ 
ten as 


d 3 p 
T \ 3/2 

— J exp|-|(n-n ') 2 }>. (44) 




I 


^(n + n , )7 1 (n — n')^ = (n^n') 2 


(39) 


Using the anisotropy vector held 7 1 , it follows that 
(n^p) 2 = Wf- (45) 

and similar to Eqn. (44) we can find the matrix elements 
of the ansatz 


Notice also that the k = 0 component of 
( Qaa — QabQbbQba ) is zero 9 so that doesn’t 
contribute to the k = 0 part of the reduced density 
matrix, which becomes 


RE = \ af + 


(46) 


so that, using a saddle-point integration and the fact that 

(7 1 ) 2 = -1, 


Pa [7^4,0 > n A, 0] 
oc 


exp [-|p||C||(n A , 0 - n' Ai o) 2 ] (40) (v\ e \ x ) ~ exp 
x exp [-|(p_LCx - p||C||) {n A pl 1 n’ A 0 ) 


-y-(* - v ) 2 - 

2 a 2 a + b 


where the effective moment of inertia is given by 
1 = \j d 2 xd 2 yQ(x,y). 


( 47 ) 

Equating this result with the reduced density matrix, 
Eqn. (40) gives the entanglement Hamiltonian 


( 41 ) 


hL = 


1 


_ _ v T 2 

tos o T „ „ / / • 


21 P I| C H 


o;/3 


{p±C± - P||C|| 


ct>/3 


84p||C||(2p|| C || - p_lc_l) 


C. The Tower of States 

We have come to the centerpiece of this paper: the 
entanglement spectrum of the k = 0 component of the 
reduced density matrix, Eqn. (40), displays the tower of 
states. 

To prove that, we need to find the entanglement Hamil¬ 
tonian H^ os such that 


(TrZ/ 7 1 ) 2 . 

( 48 ) 

Finally, we are in a position to find the low-lying states of 
the entanglement spectrum. For that we need to find the 
representation of the angular momentum algebra of the 
particle on the 3-sphere. This algebra is isomorphic to 
0(3) x 0(3), which can be made explicit by defining 24,25 


p A [n A ,o,n' AO ] = (n A , o|e Ht °’\n' Afi ). 


(42) 


The density matrix Eqn. (40) describes the motion of a 
particle on a 3-sphere. The dynamics of a particle con¬ 
straint to move on a circle should be expressed in terms 
of the angular momentum , 15 


Lap — BaPfi TlpPa 


(43) 


•Pi 2 ^ 2 ^7 ^ La% 

Ni = 2 (~2^ijkLjk Lai 


such that their commutation relations are 

Jj\ — l€,.jk dh 
[Ni,Nj] = ie ijk N k 
[JuNj] = 0 . 


(49) 

(50) 


( 51 ) 

(52) 

(53) 
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Note, however, that J 2 = N 2 = jL 2 ^, so that the rep¬ 
resentations of 0(4) can be labeled by a total angular 
momentum j. States in this representation are labelled 

by 

|{j}mi?n 2 ) (54) 

where ?n,; = —j, — j + 1,... j are the quantum numbers 
for each 0(3) subalgebra. For each value of the total 
angular momentum there are (2 j + l) 2 states. Note also 
that in this representation the anisotropy term becomes 
relatively simple, 

i(Tr Ly 1 ) 2 = (W) 2 . (55) 

The eigenvalues of the entanglement Hamiltonian Eqn. 
(48) are thus 


This corresponds to the degeneracies observed in Ref. 11 . 
For each value of total spin j, there are (2 j + l) 2 states 
in the entanglement spectrum, split into 2 j + 1 states 
with different value of m 2 . The structure for each value 
of j, which is shown in Fig. 2 and 3 of Ref. 11 . The 
fact that for integer j the highest energy state is single 
degenerate, suggests that this highest energy state has 
?n 2 = 0. Consequently, the prefactor in front of the (Vi) 2 
term is negative for the model they studied. 

We can compare our results with the square lattice 
antiferromagnet, where the order parameter is an element 
of the 2-sphere symmetry. In that case the entanglement 
spectrum has a (25 + 1) degeneracy 8 . Only in the case 
of the triangular and Kagome antiferromagnets, where 
there is an 50(3) coplanar instead of a collinear 5 2 order 
parameter, we find a (25 + l) 2 degeneracy. 


K 

,7712 


7 ^n c n 


UU +!)) 

{P±C± - P||C|| 


2/p||C||(2p||C|| - p±c±) 


(?n 2 ) 2 (56) 


The attentive reader might ask: what happened to the 
k ^ 0 components of the reduced density matrix? Here 
we can be brief, as this part has the shape of the density 
matrix of an harmonic oscillator, 


Pa[k, w'} oc exp 


— -( 7 r — — 7r') — -( 7 r + 7r , )Af 2 (7r + it') 


r 


(57) 


We can thus directly apply the results from Ref. [9] to 
find that the entanglement Hamiltonian for the spin wave 
spectrum is just that of free Goldstone bosons, the dif¬ 
ference only being in the anisotropy of the spin wave 
velocities. 

The full entanglement Hamiltonian is thus given by 
the tower of states contribution of Eqn. (48) and the 
spin-wave part, 


T 2 

h e = L ^i 

i°s 2/i 


(TrLy 1 ) 2 

2 1-2 


£ 

e,a 


/ y e a a e,ot a e 


(58) 


where we have, for compact notation, absorbed the p and 
c-dependence into the parameters I\ and / 2 , and a runs 
over the three different spin waves with energies e a . 

How do the parameters I \, / 2 and e a depend on the size 
of the subregion? Here we follow the results from Ref. [9] 
again. Consider a two-dimensional torus divided into two 
cylinders A and B by making cuts at x = 0 and x = i. 
On such a symmetry, the lowest energy eigenstate for the 
spin-wave spectrum scales at e ~ (log (,/a) , where a is 

the UV cut-off. The moment of inertia parameters Eqn. 
(41), on the other hand, scale as 

1% ~ t\og{£/a) (59) 


which makes the gap to the first excited tower state an 
order of magnitude £ smaller. Therefore, the tower of 


states entanglement spectrum is separated from the spin 
wave entanglement spectrum. 

In conclusion, the low energy entanglement spectrum 
has the same structure as the tower of states. Therefore, 
the degeneracies of the entanglement spectrum allow for 
a determination of the symmetry of the relevant order pa¬ 
rameter, similar to the idea proposed that the entangle¬ 
ment spectrum allows for a determination of topological 
order. 7 


D. Entanglement entropies 

Finally, let us mention the consequences of the 50(3) 
order parameter to the entanglement entropy. The Von 
Neumann entanglement entropy is defined as 

S E = -Tr (p A log p A ) (60) 

or can be equivalently found as the limit of Renyi en¬ 
tropies 


=-^Tr (p n A ) (61) 

n — 1 

such that limn-*.! S n = Se- 

It turns out to be quite undemanding to compute the 
Renyi entropies. Let us therefore focus first on the tower 



of states part. Starting from the reduced density matrix where / ~ flog t/a. The reduced density matrix squared 
expressed in terms of the 7r-vectors, we can get higher is 
powers of the reduced density matrix p^[ 7 r, 7 r"] by Gaus¬ 
sian integration. Recall that the relevant part of the re¬ 
duced density matrix is 


Pa[> r 0 , 7 Tq] ~ exp 


~2^{' k o ~ irb) l Cij(no 



(62) 


Pa N, K\ 


d 3 7 


v / det ~p 

d 3 Sn 

^/detp 

3/2 


exp 

exp 


- 7To)Vy(7r 0 - n' o y - ^/( n’ 0 - tt^c^ tt ' 0 - tt'o ) 3 
- ISn l Cij6n 3 + ISn i c i j( 7r 0 - ttqY - \l (^o - - KY 


2 tt\ 

21 


(det cp) exp 


--/(ttq - n'dyciji 7T 0 - 7 i' 0 'Y 


(63) 

(64) 

(65) 


Notice that we included the spin stiffness as part of the measure on the integration over tt. In the second line, we 
defined 5n = itq ~ tt' 0 . This integration technique can be continued to find the n-th power, 

/o_\ 3(«-i)/2 

P\V o.tto] ~ n ~ 3/2 {-jJ (detc/3) _(n_1)/2 exp 


KYcij^o - ' rr'oY 


( 66 ) 


The Renyi entropy associated with this result is 


Sn = g l 0 g 


I 

27r 


3 logn 
2 n — 1 


+ ^ log det cp. (67) 


The entanglement entropy of the tower of states becomes 

S E = ^ log () T ~ T ~ log det cp ( 68 ) 

3 

~ - log t + const. (69) 


The tower of states thus introduces a logarithmic term 
in the entanglement entropy. The large-size of this term 
is, to leading order in £, independent of the anisotropy 
introduced by the coplanar order. In fact, it confirms 
the picture that the logarithmic correction in symmetry 
broken systems scales as ^4/-! log t where TV — 1 is the 
number of Goldstone modes. 

The spin wave part will contribute to the entanglement 
entropy in a similar fashion to the collinear results, which 
yields an area law Se ~ t with a nonuniversal prefactor. 


j and have ( 2 j + l ) 2 degeneracy which is lifted by the 
spin-wave anisotropy. Consistent with earlier results for 
the collinear antiferromagnet, the entanglement entropy 
counts the number of Goldstone modes of the broken 
symmetry state. 

These results point towards interesting universal be¬ 
havior. Will a system that undergoes spontaneous sym¬ 
metry breaking always hide this in its finite size ground 
state entanglement? Clearly, ferromagnetic ground 
states do not exhibit this behavior as they are known to 
be exact product states. For the antiferromagnets stud¬ 
ied here and in Ref. [9], we already knew that there was 
a tower of states associated with the symmetry breaking. 
If one would be able to prove in general, that systems 
that break a symmetry in the thermodynamic limit 28 
should have a tower of states in their entanglement spec¬ 
trum and a logarithmic correction to the entanglement 
entropy, then such a result would be invaluable to search 
and identify yet unknown broken symmetry phases of any 
Hamiltonian. 


III. CONCLUSION 

Using an SO(3) nonlinear sigma model, we have shown 
that the ground state entanglement spectrum of a copla¬ 
nar antiferromagnet displays the ’tower of states’ struc¬ 
ture. States in this tower can be labelled by total spin 
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